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Abstract 



We consider m-th order linear recurrences that can be thought of as 
generalizations of the Lucas sequence. We exploit some interplay with 
matrices that again can be considered generalizations of the Fibonacci 
matrix. We introduce the definition of reflected sequence and inverted 
sequence and we establish some relationship between the coefficients of 
the Cayley-Hamilton equation for these matrices and the introduced 
sequences. 

1 Antefacts 

Let us define the m x m matrix A m as a matrix with the first column of 
all ones, as well as the first upper diagonal, while all the other elements are 
equal to zero, that is 



< 




1 for i = 1, . . . , m 
1 for i = 1, . . . , m — 1 
otherwise. 



We can write in partitioned form 



A 



m—l 





where ei and e 2 are m x 1 vectors 



ei 



Note that the inverse of A m is a matrix B m such that 
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e 2 
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_ _ 



r b 



lm 



bij 



= — 1 for i = 2, . . . , m 
= 1 for i = 2, . . . , m 
= otherwise. 



If we evaluate the determinant developping according the elements of the 
last row we see that the only non-zero summand corresponds to element in 
position (to, 1) while the corresponding minor is a m - 1 x m - 1 identity 
matrix: it follows 

iA m | = (-ir +i . 



_ m -l_ x m-2 X _ L 



(2) 



The characteristic polynomial of A m is 

g(x) =X m -X n 

Now, just to fix notation, given n numbers 

i a i}\ = {°i> a 2, a n } 
let S% (cii) denote the symmetric functions, that is 

s 'k ( a i) = a h a i2 ■ ■ ■ a ik , 1 < k < n. 

h<i2<---<ik 

The Cayley-Hamilton equation for A m is 

\m _ \m-l ... a T — n 



Using relationships among coefficients of the Cayley-Hamilton equation and 
the eigenvalues {r^, i = 1, . . . , m} we have 



Sr(r t ) = (-l) m , 



, TO. 
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Note that the maximal real root approaches 2 for m going to infinity (see 
From the characteristic polynomial we can define the m-th order recurrence 

UW=ut\ + utl + --- + utl n>m. (3) 

With m = 2 we have the Fibonacci-Lucas sequence, with m = 3 the Tri- 
bonacci sequence, with m = 4 the Tetranacci sequence, and so on. For this 
reason we might call this matrix a Polynacci matrix (Polymatrix) . 
Note that (see [§) 

n n—ii n—ii i m -2 

E E ••• E r i r 2 >< • • • x fcrjr* 1 -"-*"-* = K+i, 

U=0i2=0 i m _i=0 

where 

Ki = V n -l + 14-2 + • • • + Vn-m+l + Vn-m, 

with initial conditions 

Vb = = V 2 = ■ ■ ■ = V m -2 = Vm-1 = I- 

2 Generalized Polynacci Sequences 

Now we want to choose the initial conditions of Un in such a way that, for 
any m, there holds the Binet form 

tfM= r » + r J + ... + r n. 

In this way we obtain generalized Polynacci sequences. The term general- 
ized stems from the fact that the Tribonacci numbers (m = 3) so defined 
bears with the Tribonacci with initial conditions 0, 1, 1 (see Q) the same 
resemblance as the Lucas sequence does with the Fibonacci sequence. 
It follows 

= tr(A« ). 

Now we are going to establish the following result 

tr(A£j = tr(A^_ 1 ), k < m, (4) 
where tr(-) is the trace operator, so that 

tr(A*J=tr(A*). (5) 
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We assume k > 0, since tr(A^) = m. We need the following result, which 
can be established with tedious calculations, 

e^ei = 0, i = 0, 1, . . . , m - 2. (6) 

Using repeated matrix multiplication on the partitioned form (Equation |l|) 
of A m the asked for trace is the sum of the trace of the resulting first diagonal 
blok plus the resulting scalar which is the second diagonal block. Starting 
analyzing this scalar we see that, using Equation |6], for h < k, this term in 
is e^A^l^ei. Since h < k and k < m it follows h — 2 < m — 2 so that 
using again Equation || with m replaced by m — 1 we have e^A^I^ei = 
0, V/t. 

As for the first diagonal block the resulting expression in Aj^ is a rather 
messy sum one term of which is A^ l _ 1 . But since we are interested just in 
the trace, using the fact that for a matrix B and vectors u and z we have 
tr(Auz') = tr(z'Au), we get that the trace of this block is t^A^-J plus 
the sum of terms of the form 

* II e&Aj£?ei, (7) 
i=i 

where 




that is a = if k is even, and a = —1 if k is odd, and q are constant. Since 
the greatest power of A m _i in Equation || is k — 2 and k — 2 < m — 2 again 
invoking Equation |6| all these terms turn out to be equal to zero. From this 
we obtain Equation ||. 

It follows that we can determine the initial conditions recursively from Lucas 
initial conditions (m = 2), obtaining, for example, for generalized Tribonacci 
(A001644 in §) 3,1,3; for generalized Tetranacci (A073817 in §) 4,1,3,7; for 
generalized Pentanacci (A074048 in ||) 5,1,3,7,15; for generalized Hexanacci 
(A074584 in @) 6,1,3,7,15,31. 

But we can also obtain a general formula. As we know = m, = 1. 
We are going to prove that 

jj(m) = y - 1, i = 2, . . . , m - 1. (8) 

This will be done easily by induction on m. This is true for m = 2 since 
is L2 which is 3 = 2 2 — 1. Assume that the claim holds for m: then we 
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have to show that 



jj(m+l) =2 i_ 1 



Now, for i = 2, . . . , m — 1, 

U\ m+l) = tr(A^i) 
= tr(Ajj 
= U$m) 
= 2 i - 1, 

because of Equation f|. For i = m we have 

+ ^ = tr(A m+1 ) 
= tr(A-) 

But Un^ is the first determined by the recurrence, so 

m— 1 



t4 m) = E^ (m) 

=0 

+ l+£(2 i -l) 



i=0 

m— 1 

m ' 

8=2 

2 m - 1, 



since 

m— 1 

E (2 i - 1) = 2 m - m - 2, 

i=2 

and so the claim is proved. 

Now we are going to derive a closed form of the ordinary generating function 
(ogf). So let G(x) be the ogf of ui m) 

G(x) = clq + a\x + ci2X 2 + CL3X 3 + • • • 
where a« = U^ m \ Multiply both sides by x i , i = 1, . . . , m. Then 

G(x)(l — x — x 2 — x 3 — ... — x m ) = ao + (ai — ao)x + (012 — a\ — ao)x 2 
+(a 3 - a 2 - ai - a )x 3 H h (a m - a m _i ai - a )x m H 
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Because of the recurrence relationship all the coefficients of x J , j > m are 
equal to zero. Now insert the initial conditions a,i, i = 0, . . . , m — 1 and we 
are left with 

G(x)(l - x - x 2 - x 3 - ... - x m ) = m + (1 - m)x + (2 - m)x 2 H h 

(m-2 \ 
2" 1 " 1 - 1 - (2 4 - !) - 1 - m ) xm ~ l - 

The last summand turns out equal to —1. Then we can conclude 
m — (m — l)x — (m — 2)x 2 — • • • — x m_1 

Note that we can allow for negative subscripts, following and generalizing 
|J. If we have a general m-th order linear recurrence 

Ut ] = siUl™\ + s 2 ui™l + • • • + s m utL n>m (10) 
then we can define 



rrM _ s m-l TT (m) s m-2 rf (m) _ s l rr (m) , 1 rr (m) 

n ~ ~71 -I™- 1 ) s U -(n-2) "' s _ u _(„_ m+ i) "I" ~ ^_( n _ m )- 



fill 



3 Reflected Sequences 

The reflected polynomial ((||, p. 339]) of g(x) is 

g R ( x ) = - x m - x m ~ l - x m ~ 2 x + i. (12) 



The roots are the reciprocals of the roots of g(x), that is |fr j> an d so g 
is the characteristic polynomial of matrix B m . 



Definition 1 The reflected recurrence Ut ] of recurrence [J is the recur- 
rence with characteristic polynomial which is the reflected characteristic poly- 
nomial and with initial conditions such that the coefficients of the respective 
Binet forms are the same. 

It follows 

fj(m) _ _r/M _ fj( m ) fj(m) , fAm) > 
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and 

(Am) _J_,J_i — trfB n ~1 

' 1 '2 'rrt 

Now we are going to evaluate tr(B^), for n = 0, 1, ... m — 1 so that we get 
the required initial conditions. Of course tr(B^) = m. We know what is the 
expression for B m , so tr(B^) = —1. Now consider what happens when we 
perform matrix multiplication, starting from B m B m : when we postmultiply 
by B m then the columns of the matrix to the left are shifted to the left 
by one place and the last column is the linear combination of the columns 
of the matrix to the left with coefficients {1,— 1,— 1,...,— 1}. It is easy to 
see that for B m B m the last column is {—1, 2, 0, ... , 0}. Then all diagonal 
elements are zero except for the element in the next-todast column which 
is —1, so that the trace is —1. If we go ahead and consider B^ we see 
that the last column is {0, —1,2,0, ... ,0} and the only non zero diagonal 
element appears in column m — 2 and it is equal to —1, so the trace is —1. 
Repeating the same reasoning we come to B™ -1 : here the last column is 
{0, . . . , 0, — 1, 2, 0} and the only non zero diagonal element appear in column 
2 and it is equal to — 1 , so the trace is — 1 . So we have proved that the asked 
for initial conditions are 

u — m, u 1 — u 2 — — u m -\ — 

Of course either using the recurrence or going on with the multiplication 
process it turns out that Um^ = 2m — 1. If for example m = 4 we get 

t/ (4) = 4, &?> = -1, Of = -1, Ut ] = "I, 

while 

= 4, t4 4) = 1, = 3 , L#> = 7. 

Now we are going to derive a closed form of the ordinary generating function 
(ogf) for the reflected recurrence. So let G{x) be the ogf of Un^ 

G(x) = ao + a±x + ci2X 2 + a%x 3 + • • ■ 
where a« = u!> m \ Multiply both sides by x l , i = 1, . . . , m. Then 

G(x)(l + x + x 2 + -... + x m ~ x - x m ) = a + (ai + a )x + (a 2 + oi + a )x 2 
+(03 + a 2 + ai + a )x 3 H h (a m + a m -\ -] h ai — ao)x m H 
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Because of the recurrence relationship all the coefficients of x J , j > m are 
equal to zero. Now insert the initial conditions a,i, i = 0, . . . , m — 1 and we 
are left with 

<5(x)(l + x + x 2 + ... + x m_1 - x m ) =m+(m-l)x+(m- 2)x 2 H 

+(2m - 1 - 1 1 - m)x m ~ 1 . 

The last summand turns out equal to 1. Then we can conclude 

~. . m + (m — l)x + (m — 2)x 2 + • • • + x™ -1 
W ~ 1 + x + x 2 H h x" 1 " 1 - x m ' ^ ^ 



Note that for generalized Polynacci sequences Equation 11 becomes 



77 M _ 77M 77O) _ 77O) 1 77M 

U -n - U -{n-l) U -(n~2) "' U ~(n-m+l) U -(n-m) ' ^ 

From Equation || we get easily 

TT (m) _ TT (m) TT (m) TT {m) TT (m) 

u -l — u u l " ' u m-2 u m-l 



1. 



In the same way 



77O) _ 77M 77M rrM , 77M 

= -1, 

and in general 

f7_(i_i) = -1 i = 2, ...,m. 
From this it follows that 

^ m) =^n } - (15) 

Reflected Tribonacci is A073145, reflected Tetranacci is A074058, reflected 
Pentanacci is A074062 in 0. 



4 Inverted Sequences 

Another related sequence Ur, 
generating function G{x) as 



Another related sequence is obtained in the following way. Define its 



G(x) = -g(-). 
X \x 
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Then we get 

A . , 1 + 2x + 3x 2 + • • • + (m - l)x m ~ 2 - nix™' 1 
Glx) = ~ , . 

v ' l + x + x 2 ^ — + x m ~ x - x m 

Note that the numerator is the derivative of the denominator and also that 
the denominator is the same as in G(x). Using the Rational Expansion 
Theorem for Distinct Roots in g, p. 340] we get easily the closed form 

(Am) _ 1 1 1 



r n+l n+l n+1 
1 1 1 2 ' m 

= -U&l (16) 

The recurrence is then 

Mm) _ _f T (m) _ f T {m) _ _ fJm) f T (m) 
u n ~ u n—l u n-2 ' ' ' u n-m+l ~r u n-mi 

with 

= 1,U^ = 1, ...,U^ 1 = l-2m. 
tji 3) is sequence A075298 in §. 

5 Cayley-Hamilton Equation 

Now consider the Cayley-Hamilton equation for A 1 ^ 

A™ - ci n) A"( m ^ + • • • + (-l) l c\ n) A n Jr~ l) + ■■■ + (-l) m 4 n) I = 0, 

(n) 

where c\ is the sum of the determinants of the principal minors of order i 
of A^. Immediately we have 

c l ^ = t r (Am) = 

„(n) _ I a \n _ (_l\(m+l)n 
L m ~ I -"-ml — \ 1 ) 

Note that if m is odd c*m = 1, while if m is even = (— l) n . 
Then 

J_ + J_ + ... + J_ ^-i (C) 

r l r 2 r m lli=l r i 



J m— 1 



(_l)(m+l) 
tr(B» 

77 M 



77 
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From this it follows 

& = (-l) (m+1) ^ m) , (17) 

that is c^Lj = Un if m is odd, and c£{-i = (—l) n Un if m is even. Note 
that if the generating function of is G(x) then the generating function 
of (-l) n ui m) is G(-x). 
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